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Abstract 

Equations for the mean and the turbulent quantities for compressible turbulent 
flows are derived in this report. Both the conventional Reynolds average and the; 
mass- weighted, Favre average were employed to decompose the flow variable into a 
mean and a turbulent quantity. These equations are to be used later in developing 
second-order Reynolds stress models for high-speed compressible flows. A few recent 
advances in modeling some of the terms in the equations due to compressibility 
effects are also summarized. 
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1. Introduction 


This report describes the progress that was made during the first phase of an 
effort to develop new second-order closure models for compressible turbulence. 

Compressible turbulence modeling is an essential element of many problems 
of practical interest, such as external aerodynamic calculations, the design of en- 
gine components and noise reduction. Initially, based on Morkovin’s hypothesis 
(1964), the direct extension of incompressible models were used in the calculation 
of turbulent flows at moderate Mach numbers. This practice has enjoyed a con- 
siderable success in the calculations of wall shear layers in the past. However, it 
failed to predict adequately the reduced growth rate of high speed shear layers, in 
which the compressibility ('fleets are more prominent. As a result, assuming that 
the turbulent fluctuations undergo thermodynamic processes such as isentropic or 
isothermal processes, Oh (1974), Rubesin (1976) and Vandromme (1983), among 
others, added compressibility corrections to the baseline incompressible models and 
produced mixed levels of success. 

However, with the recent interest and need to develop high-speed flight tech- 
nologies, models that predict correctly the dynamics of complex compressible turbu- 
lent flows are needed to measure up to the advances in other areas such as CFD and 
aerothermodynamics. To this end, the present task explores second-order models 
with an explicit account of compressibility effects. With this in mind, we first derive 
the equations for the mean flow and turbulent quantities within the framework of 
second-order modeling by using two basic approaches, i.e., the Reynolds average and 
the Favre average. These equations and some related matters that were observed 
along the way are described in this report. Two recently proposed arguments about 
the influence of compressibility are also described. These include the concept of 
dilatation dissipation and the description of the changes of turbulence structures by 
flow instabilities. It is concluded that the instability wave description of turbulent 
large structures complements the second-order modeling methods and, together, 
the two approaches may be able to satisfy a broad range of needs in engineering 
calculations involving compressible turbulent flows. 

2. Basic Compressible Equations 

The equations that govern the flows of compressible fluids are the Navier-Stokes 
equations. The equations are 
Equation of conservation of mass 


P,t + ( P u i),i = 0 

Equations of conservation of momentum 

(P u i),t T = & ij,j 


( 1 ) 

(2) 
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Equation of conservation of energy: in total energy 


(,pE\t "h (pEUi) t i <Ji t i 


( 3 ) 


or 


Equation of state 
where 


{.pH u i\i P,t d" ( T ij u j\i Qi,i 
P = pR T 


(4) 

(5) 


1 Vi 

E = e + -UiUi , H — e-\ h -u.-u; , e = C v T 

A p A 

° V = ~P &ij + T ij » r ij = 2 P s ij~ P* u k,k^ij , sy = 


P o P Pbulk i ?i — h J 1 ,' 


With the expressions for the transport coefficients {p , //* , k} obtained experi- 
mentally, the above equations form a closed system with appropriate boundary and 
initial conditions. Note that according to the second law of thermodynamics, the 
transport coefficients have to be greater than zero. 


3. Equations for Compressible Turbulent Flows 


Statistically, the random fluctuation of flow properties in turbulent flows can 
be decomposed into an average value and a fluctuation. Two commonly used aver- 
a gi n g techniques are the conventional and mass- weighted averages, which are often 
referred to as the Reynolds average and the Favre average, respectively. In the 
Reynolds average, the flow property, <f>, is decomposed into a ensemble mean value, 
<f), and a fluctuation, <f>", i.e., 

P = P + p" 

Ui = Ui + u" 

p ~ p + p" 

T = T + T" 

E = E + E" 

H = H + H" 

P = P + p" , P* = V + p*" 

Cv = C v + C v , Cp = Cp + Cp ( 6 ) 
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A Favre-averaged or mass- weighted- averaged mean value, <j> is defined as 



( 7 ) 


A mixed average that uses the Reynolds average for the density, pressure, transport 
coefficients and specific heats and Favre average for the other quantities is often ap- 
plied in dealing with compressible turbulence. Flow properties can then be written 
as 

P = P + p" 

Ui = Ui + u'i 

p = p + p" 

T = T + T' 

E = E + E' 

H — H + H' 

P = fl + p" , p* = p* + p*" 

c v = c~ v + c” , c p = c~ p + c; ( 8 ) 

Note that both decompositions are defined in term of the ensemble average and 

< = 0 , < ^ 0 ( 9 ) 

= 0 , ± 0 (10) 

This is one of the most important differences between the Reynolds and the Favre 
averages. The relations between these two averages are 


<t> — 4> = 4> ~ 4>" — <t>' = - p"4>' / p = - p"<t>" / p (11) 


3.1 Mean Flow Equations 

3.1.1 Reynolds-averaged Equations 

The Reynolds- averaged mean equation for the conservation of mass can be 
obtained by taking the Reynolds average of the instantaneous equation of the con- 
servation of mass. The mean equations for the conservation of momentum and 
energy can be obtained by the same procedure. 

Mean Continuity 

P,t + (pui + p"u'?) }l = 0 (12) 
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Mean Momentum 


P U >,t + (p"u' l '),t+ (p Uj + p"u'l) Uij + (pu”u” + p"u{ , u'J + u j p"u'')' j 


— { (P + P* u k,k + P* u 'k t k)&ij + 2(p s ij + p" s ij)}, 


( 13 ) 


Mean Energy 


(pE + P "E"\ t + ( pEu i + pE"u '/ + u iP "E" + E p”u” + p"E"u '! ) (i 


{ +P Ui +P"U 1 + p*U ky k Ui + P* u 'k y k U 'i + U iP*" U k,k + U k,k P*" It” 


+ P*" u k,k u i) + 2 (p s ijUj + ps” jU ” + Ujp"s'!j + l i3 p"u” + /d'u"s") }, t (14) 


or 


{pH + p"H" - p) it + ( pHui + pH"u” + u tP "H" + H p"u" + p"H"u") ti 


= {-(Qi+P* u k ,i tw, + p* u” k u'! + Uip*"u” k + u k , k p*"u'! + p*" u'^ k u”) 

+ 2 (p SijUj + ps” 3 u” + + s tj ]PWj + p"sVjU'!) },* (15) 

where 


and 


Equation of State 


E — C v T + C"T" + - (uiUi + u'-u”) 


H = c v T + c;T" + - (u t u, + «) 


p = pRT + Rp"T" (16) 

These mean equations may be used to solve for the mean quantities: 

{? i «i i P , T} (17) 

The set of mean flow equations, however, are not closed. Quantities that need to 
be known to solve these equations are 


P" u " , w"w" . P"u>j , P*"ulk , , />"£" , E”u” , , p"u” , 






U *,* U < ’ ’ P*"«M » p*" u lk u i » 

P 77 ^ , , JF , p , C~ v (18) 


"«"s? 
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3.1.2 Favre-averaged Equations 


The equations governing the Favre-averaged mean variables can be obtained 
by th e sam e process that was used to derive the Reynolds- averaged mean equations. 
Since po' = 0, the density fluctuation does not appear explicitly in the resulting 
equations. The equations are listed below. 

Mean Continuity 

P,t + = 0 (19) 

Mean Momentum 

C m),t + (puiUj + pUiUj),j = 

{-(p + ]Fu k>k + p*"u' kk + T^ k )Sij + 2 (psij + + P (20) 

Mean Energy 


(pE),t + ( pEui + pE'u'j)^ = {-(^ + pui + pu't + p"u\ 


+ P*Uk, k %i + P*u' k k u' i + u i p*”u' k>k + u kjk p*"u\ + p*"u' i u , Kk + p*u l u' k k + p* Uk'ku'i) 

+ 2 (jZ Sijuj + -ps'ijU'j + Ujp"^ + + J£ujJ- } + psijy^)},i (21) 

or 

{pH - p) (t + ( pHui + ■),» 

- i~ (Qi+]Su ktk Z i +JFu' ktk u[ + iZip*''u , k ^ + u k p t p*" u' + + JFui^) 

+ 2 {psijuj + Jls'ijU'j + + Sijp"^ + p ir u i ^ j + p u-jlCj + } ;l (22) 

where 


E = 

C V T 

+ C"T' 

+ kfcw 

r +Q' 

P"T') + 

\ (“••«< + u ', u 'i) 


and 



1 - 



. 


H = 

C~ P T 

+ C"T' 

+ =(TC”pP 
P 

r +c" 

P"T>) + 

2 («i“» + u ’i u 'i) 


Equation of State 










P = 

pRT 



(23) 

These equations 

may 

be used 

to solve for 

the mean quantities: 





{p,Ui 

i P , 

T} 


(24) 


Not surprisingly, the set of mean flow equations are not closed. Quantities that 
need to be known to solve these equations are 



It was shown in this section that, in order to solve the mean flow equations, 
the turbulent moments of order up to three must somehow be modeled. The main 
objective of this study is to develop rational second-order modeling procedures 
for compressible turbulence. The transport equations for the second-order tur- 
bulent moments are derived in the following sections. Higher order moments are 
represented as functionals of the mean flow and the second-order moments in the 
second-order models that will be developed later in the second phase of this project. 


3.2 Equations for Turbulent Fluctuations 
3.2.1 Reynolds-averaged Equations 

The equations for the turbulent fluctuations are needed in the derivation of the 
equations for moments of various order. The fluctuating parts of the equations of 
conservation of mass, momentum and energy were obtained by subtracting the mean 
equations from their corresponding instantaneous equations. The same procedure 
is also applicable to the equation of state and the thermodynamic relations. 
Fluctuating Continuity 


p"t + (pu '■ + p"ui + p'V/ - = 0 (26) 

Fluctuating Momentum 

(A:^Af+Ai+K),« + 

+ Uj(p"uV - ~p^i) + A"u" - TTJyj + P Tiiu'j + P Uju'! + u,Ujp"}j 
= {~(p" + P* u k,k ~ P*' u k,k + P* u k,k + P* Uk,k)6ij 

+ 2 (p" s ij - p" s ij + P s "j + p"sij},j (27) 

Fluctuating Energy 

(pE" + Ep" + p"E" - ~p*W), t + {p{E"u'l - E"u") + Ui(p"E" - ~p ir E n ) 

+ E{p"u'l - + p"E"u'l - /X + p TiiE" + -pEu'l + UiEp"} ti 

= ( - U"+K 

+ u kt t(n‘ u" - /J *" u") + #** «*,*«"- + ff' uiUt.k] 

+ 2 Wltf - Wj) + - T^Tj) + - A") 
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+M u j s ij ~ + l L u j s ij + !“> ijUj + 


( 28 ) 


or 


(pH + Hp + pH — p"H") r t + {p(H"u'l-H"u'l) + u l (p"H n -~pHH) 


+ H(p"u'l - p"u’!) + p"H"u'; - p"H"u'{ + p u t H" + pHu'l + u.Hp"}^ 


+ u k , k (p* u { -p*"u'?) + p* u , l iik u'!~p*' , u , l tk u , ! + p*u i ul k +p*u kjk u , l + p* UiU*, fc ] 


+ 2 Wlrf - sVtf) + u^'s'lj - p"s'r) + s^p"^ - p"u'j) + p'u'jsij 


- p"u'js" j + p Ujs'-j + psiju" + u j s ij p"]} i i 
Fluctuating Equation of State 


p" = R(pT" +p"T + p"T" -p"T") 
Fluctuating Thermodynamic Relations 


e " = c”t + c v t" + CX - C':T" + u,u'; + ^'v/ - <<) 

i 


h = c;t + c p t" + c;t" - c”t» + ay? + ±(«;v/ - ^y>) 

3.2.2 Favre-averaged Equations 

Fluctuating Continuity 


( 29 ) 

( 30 ) 

( 31 ) 

( 32 ) 


f>"t + {p"«i + (p +/»")«'},, = 0 


( 33 ) 


Fluctuating Momentum 


{pu'i + ^'(ui + u'i)} >t + {p"u l Uj + (p + p")u' i u' j - pu' i u' j + (p+p")u i u' j + (p + p")u' i Uj} J 

= {-\p" + H*"u'k,k - P*"u' k ,k + P*( U 'k,k - u)^k) + P*"u k ,k] 6 ij 

2 [p Sjj — p" s- + p( s ij — s'ij ) + p ■ s i>]},j (34) 

Fluctuation Energy 

{p"E+(p+p")E'} it + {p"Eui+(p+p")E'u' i -pE'u' i + (p+p")u i E l + (p+p")u , i E} i i 

= {-\p( u 'i - u' i )+p"u i +p"u' i - p"u' i +]F(u' kik u' i - + Ui(p*"u' k k - p*"u'l k ) 
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+ u M (p*"u' - p*'X) + H*"uy kk - + p*Ui(u' M - u' M ) + /i*wjfc,jfc(u- - w') 

+/i* «»«*,*] + 2[p(s^u' - s'ijUj) + Ujiy's'ij - p"s-j) + Sij{v" u 'j ~ P" u 'j) + fi"u' J s' ij 

4 ^j S ij ”b 4 u j{$ij ^ ij ) 4" f- ISij^Uj 4" u j$ijp ]},i (35) 

Fluctuating Equation of State 

p" = R [p"T + (p + p")T'] (36) 

Fluctuating Thermodynamic 

s’ = Cf+ar+cr-^-hf^+cjprj+Miu'+hu'u;-!^;.) (37) 

P ^ 

h' = c;'r+^T'+c;'T'-^-i(T^+^ (38) 


3.3 Turbulent Moment Equations 
3.3.1 Reynolds-averaged Equations 

Since 4> n = 0, the first-order moments of all the Reynolds-averaged 

t urbul enc e fluc tu ation s are zero. Equations for the second- order moments: 
P u i u j'> pE (t u p ff u" and turbulent kinetic energy, are derived below. 

An equation for the total momentum exchanges, pUiUj can be derived by multi- 
plying the momentum equations in m by uj, exchanging the indices and adding the 
resulting equations together. An equation for ~p UjUj can be obtained in a similar 
way. 

(p UiUj) it + Ui(p"u'j) it + Uj(p"u") it + (p UiUjU m ) t m + (uiUjp"u'J n ) >m + U i(pu" U 1 ^) >m 

+ «>(?<<), m 4«y(« m Af), ra + 


= -Uj P>i - u iPj - Ui(p* u kjk + p*"ul k )j - Uj(p* u ktk + ll*"ul k \i 

+ 2 Uj(pSj m + p Sj m ) m 4 ' 2 Uj(pSi m + p"S{ m \m 


(39) 


An equation for p u"u ", the turbulent Reynolds stresses, may be obtained by sub- 
tracting the equation for p u t u 3 from the one for pu t Uj. The resulting equation 
is 


9 



(pu"u" + p"u' i 'u'^) t t + p"u'lui,t + p"u"u jt t 


+ Ui^pu'ju'^ + u m p"u'j + p"u , ju'J n ) + + U m p"u '! + p"u , !u'^ n ) 


+ 0 u m u'/u'j + pu'fu'fu'Jn + u m p"u'lu'l + p"u"u"u " ) >T 


= - u "p”i - u iP"j - u 'j(p* u ik + P*"u k ,k + P*"u'l tk ),i 


U i (P* U k,k + P* U k,k + p* u 'k,k)j 


+ 2 K(Km + P"*im + p"*L), m + + p"s jm + p” S H jm \ n } (40) 

The continuity equation and the following relations have been used in the process 

a (bc) iTn + b (ac) >m = ( abc) jTn + b a c m (41) 

a (bcd) tTn + b ( acd) >m = ( abcd) jm + b a (cd) >m (42) 

l 

An equation for the tu rbulent kinetic energy k can be obtained by letting i = j 
in the equation for pu'-u ", 


(pk + p"k), t + p"u'!ui tt 


+ (P u i u 'j + Ujp"u" + p"u"u")ui t j + ( p ujk + pu'jk + ujp"k + p"u"k)j 


= -<p" - u , ?(p*u% k + p*"u k>k + p*"u'l k )' i 


+ 2 <(ps" + p"s tj + p"s'l } ) tJ 


(43) 


p"u'' 


An equation for p"u ” can be derived by manipulating the equations for u" and 
p" ■ An equation for p" can be obtained from the continuity equation. It is 

P"t + u jP"j + u 'jP,j + P u j,j + P"uj,j = 0 (44) 

An equation for u" can be obtained from the momentum equations. Note that 


1 i — 1 ( 

P P + p" p 1 + c 


2 


V) = ^(l-^ + (^) 2 +-.) * =-^ + 0((^-)) 

~ P P P P p P 

(45) 
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This linearized form for 4 is used here to simplify the equation. With this simplifi- 
cation, the equation for u" becomes 



u 'i,t + u i u h + WjUij 


~(P + P u k,k + p"uk,k 4- p + 2 (ps" + n"sij + fi"s" '■)} ,j 


— =2 i~(P + P" + P* U k,k + P* Uk t k + P* u k,k + P* u k,k)^ij 

+ 2(p Sij + + p.s" + p"s'i'j)},j (46) 

The equation for p"u" can be obtained by multiplying the equation for p" by «'/ 
and the equation for u" by p" and summing up the two. The resulting equation is 


+ Uj{p"u ;•% + p"u" 3 u hJ + p"u''u h} + (p"u'; u 'j) tj + u'yjpj + yy^ 


= - y {(p" + P^"* + p*"u k ,k + P*"K,k) 6 ij 
-2 (ps"- + p"s i:j + p"s”j)}j 

- (y ) 2 { (p + p" + p*u k ,i k + P^" k + p*"u k ,k + P*"u ktk )h} 

+ 2 (p sij + ry + + p"^') jj + o((^) 3 ) (47) 

The linear approximation was used in Shih et al. (1987) in modeling a variable- 
density mixing layer and may be applicable in the modeling of compressible turbu- 
lence as well. 



One form of the transport equation for pE"u '/ can be obtained by tracking 
the following steps: (1) Derive the pEu, equation by multiplying the momentum 
equations by E and the energy equation by u t and summing up the two resulting 
equations (2) Derive the pE u; equations by the same procedure (3) Subtract the 
pE U{ equation from the pEv , equation. The resulting equations are 


(pE"u” + p"E"u'?), t + p"u'l E, t + p"E"u i>t 


+ ( p UjE"u" + ’pE"u"u" + Uj P "E"u'! + p"E"u"u" )j 


+ (pu'fu" + Uip"u'j + Ujp"u " + p"u'{u 'j ) Ej 
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= - E"(p" + /i*u" k + p*"u k>k + p*"u'l k )j 

- «X(p" + Wk,k + A*"«M + 


- «"«"(p + p" + p*u m + * + p*"T/ M + fi*"u'l k )j 

- Ujju”(p" + /i*u2 >Jfc 4- P*"u ktk + p*”u'l k ) 


- u j u j,j(p + p" + P*Uk,k + P*u'k, k + P*"uk,k + P*" u k t k) 


+ 2 [u''(ps ijU 'l + ps'; jUi + ps'lju'l + fl"s lJ U l + p"s ijU " + p"u t s' l ' J + p"u'fs^) tj 


+ E" (M'j + p"sij + p"s"j),j] - q'jju'i 


(48) 


3.3*2 Favre- averaged Equations 


In order to close the _Favre- averaged mean equations, the quantity u[ has to be 
known. An equation for u\ can be obtained by taking the Reynolds average of the 
equation for u\. The Favre-averaged mean momentum equations can be written as, 
by using the Favre-averaged mean continuity equation, 

(ui),t + UjHij = - {-{p + p*u kyk + p*"u' kk + p*u' k k )6ij 


+ 2 (psij + + p s-j) - puy^j (49) 

An equation for u[ can be obtained by ^subtracting this equation from the total 
momentum equations. The equation for u' becomes 


U'i,t + Uju'iJ + u'jUij + u'jU'ij 


~ -{-(p + p" + P*u k ,k + P*u' kk + P*"u k>k + p*"u' ktk )Sij} 


+ -{psij + ps'ij + p"Sij + p"s' ij ) tJ 


= ={Hp + P*Z k , k + P*"u' kk + p'u'^Sij + 2 (ps ij +p''s^+Jis , ij )-u , i u , j }j (50) 
Since, 


u' t = - p"u'J / p (51) 

an equation for u\ can be converted to an equation for the mass flux p"u". 
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Similarly, an equation for E' can be obtained from the mean continuity, mean 
energy and the total energy equations. The resulting equation is 


E', t + UjE'j + u'jEj + u'jE'j 


= ~~{ 9 j + q" + ip + p" + p*u k ,k + p*u' kk + 


u 


k,k 




2 __ _ 1 

+ + [is'ij + + n"s , ij ) yj + - { qj + ( pui + pu\ + p"u\ 

+ P*u k ,kUi + p*u kyk u J + p*u' kk Ui + p*i^ fc u'- 
+ P*”u k ,i fc«- + p*"u' kk ui + p*" u , kk u' i ) 8 ij 
- 2 (Wiju'i + WiX + P&ifii + WijU'i + p"s ll u' l + fl"s' tJ U l + n"s' l] u' l ) 

- pE*u'j},j ( 52 ) 


Note that, unlike ttj, E' does not appear explicitly in the Favre- averaged mean 
energy equation. However, this equation can be converted easily to an equation for 
p"E" which appears in the Reynolds- averaged mean equations. 

The equations for the Favre- averaged second-order moments can be obtained 
through the same process used in the derivation of the corresponding Reynolds-a 
averaged second-order moments. Therefore, only the final form of the equations 
will be given. 

puVi 


k 


(P u 'i u 'j),t + (pUmU'iUj + P u ' t u'jU' m ),m + p(u' t u' m U j >m + U , j U , m U it m ) 


= -[(p + P ") U ;],i - [(p + p"K],i + 2(p + p")s' ij 

- u'j(p*u kik + p*u' k k + p*"u kyk -1- p*"u' k k ) ti 

- u\(Jj*U kyk + JFu' k k + p*"u k<k + p*"u' k k ) tj 
+ 2 { u 'j(pSim + ps' im + p"s im + p”s' im ) yTn 

+ + ps'j m + p"Sj m + ( 53 ) 


(pk),t + (p Ujk)j + pu'iU'jUij + (pu'jk)^ 

= - [(p + p")^!],.' + (p + p"X,i 
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+ 2 u-{fiSij H~ /i-s- + l 1 ** s ij "h s\j) ,j 


(54) 



( pE'u 


i),< + (p UiE'u'i + pE'uy^j + puyjEj - pE'u'jUij 

- E'{(P + p" + P*Uk,k + P*u'k,k + P*"uk,k + P*"u k ,k)6ij 
+ 2 (psij 4- P 3 ij + p" s ij 4" p"s tj )},j 
U}J U 'i(P + p" + J^u k ,k + J^u'kj + p*"u k ,k + p*"u ki k ) 

Uju\(p + p" + p*u k< k + P*u' kk + p*"u k ,k +p*"u k ,k)j 
u j,j u t(p 4 ~ p" 4 " P* u k,k + P* u 'k,k + p* u k,k 4 - P* Uk t k) 


- u'iU'^p + p" + p*u k ,k + P*u'k,k 4- P*"u k ,k 4 -p*"u k ,k),j 

+ 2 v J 0s^u i + Jfs l3 u\ + ps' l] u l + Jls'ijU'i + p"s l] u l + p"s tJ u' t + p"s' lJ u l + p"s' l] u' l ) tJ 

- U,< - vlFi (55) 


Equations for the third-order moments can be similarly derived, if need arises, and 
may provide some insights into the physical mechanisms involved. In the mean 
time, however, we derive moment equations only up to the second order. 


3.4 Simplified Equations 


Equations for the mean flow and the second-order moments in Reynolds or 
Favre decompositions derived above are rather complicated compared to their in- 
compressible counterparts. These equations can be greatly simplified by assuming 
that the correlations of the transport coefficient fluctuations, p" and p *" , and ther- 
modynamic coefficients, C and C ” , with other turbulent fluctuations are negligible. 
That is, 


p4> 


< i , 


p*4> 


< i , 


Cjf 

C v 4> 


< 1 , 


TTJ 


< l 


(56) 


The resulting simplified equations for the mean flow and turbulence moments are 
given below. 


3.4.1 Reynolds-averaged Equations 

The Reynolds- averaged mean flow equations, the turbulent kinetic energy equa- 
tion and the Reynolds-stress equations become 
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P,t + (P u i + P" u "),i = 0 


( 57 ) 


Meai^Contiimit^ 


Mean Momentum 


( P u i,t + P"u"),t + ( P U i U j + UiP"Uj + PU'-U "■ + p"u"u'j + Ujp"u")j 


— { P $ij + T «j},j 


( 58 ) 


Meai^Ener^ 


(pE + p"E") <t + (pEui + ~pE"u” + Uip"E" + E p"u" + p"E"u”) ti 

= {-(li + P Ui +p"u'?) + TijUj + T ij u 'j },< ( 59 ) 

where 

E = C v T + i (uiili + u"u"), 

and 

T ij 2 fl Sij — fl* UkfrSij , Tjj = 2 /i 5- — /i* 

Equation of State 

p = pi?T + i?p"T" ( 60 ) 

pu'/u'- 


(P u i u j) f t + (P + P ,ru j^i,t + P ttu i u j,i 


— Pg + Tg mm + Ilg + eg ( 61 ) 

where 

P g = “ { (pu'jVJn + u m p"u'j + p"u'?u" ) u i>m 


+ (ptt'/u" 

+ U m p"u" + 

p" u y'm ) u i,rn } 


{ P u " u ” u m 

+ u m p"u"u" 

+ p"uyju\ ; + 

+ p"u"6j m 

- r? m «; - 

~ } 




ng = 2 p"s'' 


u 


= - r?' u" 

im u j,m 


u ft 

j m i, m. 
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k 


(pk) tt + (pujk)j + (p"k), t + p"u'l u M = P R + Tfj + I1 R + e R (62) 

where 

P R = ( - pu>" - u jP "u •! - p"uy> ) u ti] 

Tf = - ( yfk + u~yk + 7yk + pys tj - Vyj ) 

n R = 


The mean flow is described by the following quantities : 


\.P » u i ) P > T} 

Quantities that need to be known to solve these equations are 


u k,k u i » u " s "j > P"T" , P* , p , C v 


(63) 


p"u'f , u'lu'l , p'<u'>u'l , p"£" , £"u" , p"E"u'l , p"u" , 


(64) 


3.4.2 Favre-averaged Equations 

The simplified equations for Favre-averaged quantities are 
Mean^Continuit^ 

P,t + (pZi),i = 0 

Mean Momentum 


(pui) >t + (puiUj + pu'yj)^ = {-pSij + Tij + rb}j 


MeanJSnergy 


(pE),i + ( pEtii + pE'u[ ) tl 


(65) 


( 66 ) 


= {-(5. + + P< + P"“i) + fijUj + tPuj + fiju't + tPu'j },j (67) 

where 

e = cjr + - ( uiUi + u'y) 


and 


Tij — 2 p Sij p* Uk,k t>i j, r' tj — 2 p s'ij — p* u'k k Sij 
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Equation of State 


p = p R T 


( 68 ) 


pu'u' 


(P u 'i u 'j),t + 

~^ij "h “I” ^ij P,j^ i P j d~ ' r im,m u j "i" ’ r jm,m u i (69) 


where 


P F 


rp r 

A ijm 


P ( u i u m u j,m + u j u ' m u i,m ) 


~ { P u 'i u j u m + P"u'jSi m + p"u , l 8 jm - r/ m u' - rj m u' } 


nf, = 


•s 


2 p"*;, 


— T • U f - 
im j,m 


7”( ^ I, 

jm i,m 


C pk),t + (pujk)j = P F + Tfj + n F + e F - p t X + Tijju'i (70) 


where 


Since, 


P h = — p t/Ju' Uij 
T F = - { pu'jk +p"u' j - tIju'i } 


n F = p"u[ ■ 


»,* 


6 = - <J<3 




ng = ng and n F = n 


R 


(11) 

(71) 

i.e. the pressure-strain terms in the Reynolds and Favre averages axe identical. 
Moreover, if the viscous stress tensors are expressed in the form of Reynolds- 
averaged rates of deformation tensors, the following expressions can be derived 




eg and e F — e 


R 


(72) 
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These equations may be used to solve for the mean quantities: 


{ p , Ui , p , f } 


The unknowns are 


p"u\ 


u'u' , E>u\ 



( 73 ) 


(74) 


4. Some Observations 

The equations governing the mean motions of compressible turbulent flows 
have been derived. The mean motions were defined in two ways. Firstly, the mean 
is taken as the ensemble- averaged mean of instantaneous quantities. The mean 
motions are thus described by the mean density, velocity, temperature, etc. This 
is also called Reynolds average or Reynolds decomposition. The second averaging 
method, Favre-averaging, is a mass-weighted average. The mean flow is described 
by the mass- weighted quantities. Regardless of the decomposition that was used, 
quite a few correlations or moments of turbulent fluctuations have to be known to 
close the equations governing the mean flow. 

Due to the nonlinear terms in the Navier-Stokes equations, such as the con- 
vective terms in the momentum equations and the surface-force work terms in the 
energy equation, an equation for an n^-order moment always contains correlations 
of order higher than n. All these terms in the turbulent moment equations have to 
be modeled so that the mean flow and the turbulence equations form a closed system 
with appropriate boundary and/or initial conditions. Second-order modeling efforts 
involve solutions of the mean flow and the second order correlations that appear 
in the mean equations through rational modeling techniques. In the following, the 
equations for compressible turbulent flows obtained by using the Reynolds average 
and the Favre average will be examined. 

It should be noted that the form of the Reynolds-averaged turbulent mean flow 
equations is not the same as that of the compressible Navier-Stokes equations. The 
equations for the turbulent mean flow contain additional turbulent correlations due 
to the turbulent density fluctuation such as p"<f>" and p"u"(j)" and the turbulent 
fluctuation of transport coefficients such as p*"u' k k and p.*" S"-. As was mentioned 
earlier, the correlations of the fluctuating transport coefficients may be negligible 
in most problems. This, however, is not necessaxily the case for the correlations 
of the density fluctuation. Density fluctuation terms such as p"u" and p"E"u " in 
the Reynolds-averaged mean equations may be important and have to be modeled, 
Varma et al. (1974). Since these moments of density fluctuation do not appear in 
the equations for incompressible turbulence, there axe no incompressible baseline 
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models that modeling of these terms can be built upon. Therefore, the Reynolds- 
averaged equations are seldom used in the calculations of compressible turbulent 
flows. 

Two procedures that can be used to obtain transport equations for the density 
correlation terms were proposed in the previous sections. These equations may be 
used to develop modeled transport equations for these moments of density fluctu- 
ation. Note that the average turbulent mass flux p"u" are different from the mass 
flux fluctuation, (pui )" , since 

(pu,)" = pu " + p"ui + p'u'l - p"u" (75) 

On the other hand, Favre averaging eliminates moments of the density fluc- 
tuation and the resulting equations bear much more resemblance to those in in- 
compressible flows. This is attributable to the use of mass-weighted quantities : d>, 
since 

pif}' = 0 where <j>' = <j> — </> (11) 

In light of this analogy and Morkovin’s hypothesis that the turbulence structure is 
unaffected by compressibility as long as the fluctuation Mach number is much less 
than unity, there was common optimism for the extension of incompressible models 
to high Mach number flows. This practice has enjoyed a considerable success in 
the calculations of wall shear layers in the past. The success, however, was not 
shared by the predictions of high speed free shear layers, even when the conditions 
satisfied Morkovin’s hypothesis. The extension would also fail in the presence of 
shock waves and expansion fans where the dilatation is large. One prominent fea- 
ture that extension models fail to predict is the reduction of the growth rate of 
compressible free mixing layers as the Mach number increases. Oh (1974) consid- 
ered the pressure dilatation terms in the turbulent kinetic energy equation. These 
terms were modeled using an eddy shock wave hypothesis. Oh argued that the 
convecting turbulence structures in free shear flows provide the key mechanism to 
produce pressure-dilatation and predicted the reduced growth rate of high speed 
mixing layers. 


The Favre-averaged mean equations can be reduced to a form similar to that 
of the Navier-Stokes equations in Favre-averaged mean quantities. With an eddy- 
viscosity type of model, an operational compressible Navier-Stokes solver can then 
be extended to include “turbulence ” by simply (1) neglecting fc ands(~ and (2) 

replacing the molecular viscosity by an equivalent turbulent eddy viscosity without 
modifying the main structure of the codes. This is generally the case in the imple- 
mentation of eddy- viscosity turbulence models to existing CFD codes. It should be 
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remembered, however, that the solutions thus obtained are mass-averaged quanti- 
ties, rather than ensemble-averaged ones, since 


<p — <p — (f) 1 


p"<j>" 

P 


( 11 ) 


pu^u'j — p u"uj 


p"u"u'l 


p"u" 


p"u" 


(76) 


The differences between Reynolds- and Favre-averaged turbulent quantities may be 
less than 10% for low supersonic flows. Nevertheless, it may not be negligible for 
flows of higher speed, say for M > 5. Therefore, it is essential to be consistent not 
only in the equations that are used but also in the comparison of results. For in- 
stance, Favre-averaged turbulence equations must be used with the Favre-averaged 
mean equations and the solutions be converted to Reynolds-averaged quantities to 
compare with ensemble averaged experimental data. 


Note that <f>' ^ 0 and the Favre- decomposition looses its advantage over the 
Reynolds- decomposition in cases in which density does not appear, for instance, 
the surface force terms in the mean momentum and energy equations. There are 
a variety of ways to handle these moments. Rubesin (1990), Vandromme et al. 
(1983), Speziale et al. (1991) , among others, simply neglected S'- and u' k k terms 
in the mean momentum equations. That is, 

T{j 2 pSjj p* u k k 8{j ( 77 ) 

Sarkar and Balakrishnan (1990), instead, use Reynolds decomposed velocity in the 
constitutive relation 

Uj = ZpSij - p* u k ,k^ij (78) 

In this case, a model is needed to convert Favre- decomposed velocity to Reynolds- 
decomposed velocity, since 


p"u'l 

«*-«< = - ~r~ (79) 

p 

The surface force terms in the mean ener gy equation can be similarly simplified. 
Rubesin (1990) argued that the S' t] , u' kk and terms in the equation can sim- 
ply be omitted as in the mean momentum equations. Sarkar and Balakrishnan 
(1990) used Reynolds-decomposition in all the surface-force work terms. Speziale 
and Sarkar (1991) used Reynolds- averaged velocity in the shear stress and Favre- 
averaged quantities otherwise. All but Rubesin’s (1990) proposals used a model for 
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u\. Sarkar and Balakrishnan (1990) and Speziale and Sarkar (1990) used a gradient 
transport expansions, i.e., 


P 


P,i 


(80) 


A more refined analysis can be found in Taulbee (1991), in which a modeled trans- 
port equation was used for the mass flux, uj. 

Note that there has been a lack of consensus in the identification of the com- 
pressibility effects on turbulence, quite aside from the modeling of compressible 
turbulent flows. For example, Nichols (1990) and Grasso and Speziale (1989) calcu- 
lated separated flows and considered the effects of compressibility in the turbulent 
kinetic energy equation. They chose completely different terms in the equation 
to represent the dominant effects of compressibility. Nichols (1990) argued that 
the “turbulent velocity-density dissipation” terms, 2 v.ip"u" Sij , are the most sig- 
nificant terms; Grasso and Speziale (1989), on the other hand, proposed that the 

du" 

pressure-dilatation terms, p" and the product of the pressure- gradient and mass 

fluctuations, — -r 1 - play a strong role in high speed flows with adverse pressure 

gradients and could not be disregarded. Both Grasso and Speziale (1989) and 
Nichols (1990) used a ID isentropic approximation to model the relations between 
the density and velocity fluctuations. Nichols (1990) applied his model to a vari- 
ety of flows, which included a supersonic compression ramp flow similar to that in 
Grasso and Speziale (1989). For the compressible ramp flow, both models obtained 
some degree of sucess in the prediction of the skin friction coefficients and it is 
difficult to identify the dominant compressibility effects. 


5. Recent Advances 


Recently a rather new concept called dilatational dissipation was proposed by 
Sarkar et al. (1990) and Zeman (1990). For p bu i k = 0 and fi = constant , the 
dissipation rate 


P e 


t-'V.' 


(81) 


can be represented by 


pe = H( 


u'-'u" 

l I 


+ j K.tY ) = *. + <« 


(82) 


where is the vorticity vector and e s and e c denote the solenoidal and the dilata- 
tional part of the turbulent dissipation. The newly identified dilatational dissipation 
term accounts for the viscous dissipation of turbulent kinetic energy due to volume 
fluctuations. Both Sarkar et al. (1990) and Zeman (1990) argued that the solenoidal 
part of the dissipation assumed its standard form as in incompressible cases and 
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that the dilatational dissipation was proportional to the square of the Mach number. 
Zeman’s (1990) analysis was based on the existence of turbulent eddy shocklets in 
high speed mixing layers. On the other hand, Sarkar et al. (1990) used an asymp- 
totic analysis. Model predictions from both work were very encouraging in the 
prediction of the reduced growth rates of compressible mixing layers as a function 
of the convective Mach number. Dilatational dissipation appears to be among the 
direct consequency of compressibility effects and has to be properly accounted for 
in modeling compressible turbulence. 

Sarkar et al. (1990) indicated that e 9 is rather insensitive to the change of 
compressible indicators in moderate Mach number turbulence and modeled e a using 
a traditional incompressible form. For other cases such as for flows of higher Mach 
number a refined model may be needed. An equation of e a (= Jlu/fuf) can be 
obtained by subtracting the equation for uquq from the cCju,’, equation. The resulting 
equation is 


€s,< + U j € »,j — C .i U j,j — 2[i WiWjUjj 


2 •pu'yj&ij + 2/1 uju"u" j + 2 JIu'/u'Juij - 2e s Ujj + 2 JIu'/u'ju'lj 
i 


where 


+ ^ r " T " 

1 p *j kPiPj 


p'r'L 1 

• - — r" n n T" 

-2 T kp,PJ -2 r kjP T kp,pj 


+ ( = )b r kj T kp,p ~ r fcj(= T k P ,P ~ r k } (^l),J T kp t p ) 


r k] = ^ ( U k ,j ~ My*) 


(83) 


and S denotes entropy. The approximation (45) was also invoked. The equation 
suggests that the dynamics of c s in compressible turbulence may be very different 
from that of the turbulent dissipation rate in incompressible turbulence. The equa- 
tion also indicates that in compressible turbulence, the effects of the thermodynamic 
states of the flow system may be important to the evolution of e 3 . 

Another school of though on compressibility effects focuses on the changes of 
turbulence structures at high Mach numbers. It was proposed that turbulent ener- 
getics evolve with preferred modes of perturbations as the Mach number, or com- 
pressibility effects, increases. These energetic turbulent structures were interpreted 
as physical manifestation of evolving instability. The selective amplifying effects of 
compressibility were mentioned in Sandham and Reynolds (1987), Ragab and Wu 
(1989), Lele (1989) and Morkovin (1990). Due to the communicability problem be- 
tween interacting elements in supersonic flows, structures that are highly efficient in 
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extracting energy from the mean flow at low Mach numbers no longer prevail as the 
Mach number increases. They are replaced by structures associated with the modes 
generated by compressibility effects. Therefore, the quasi-2D vortical structures in 
low speed mixing layers may be more susceptible to the communicability problem 
than the highly 3D structures in bounded shear flows. Linear stability analysis 
also shows that for even higher M c , the amplified modes travel supersonically and 
obliquely and thus may generate only limited mixing. Morris et al (1990) applied 
linear analysis and predicted quantitatively the reduced growth of supersonic free 
mixing layers without any adjustment of operating conditions. Liou (1990) devel- 
oped three closure models to predict the evolution of incompressible mixing layers 
at the large scale using a weakley nonlinear theory. The models predicted the mean 
flow field as well as the unsteady large-scale turbulent motions. This approach 
is rather unique and appears to be very promising in predicting the unsteadiness 
of turbulent structures. Note that a reasonable modeling of the unsteadiness and 
spatial intermittence of turbulent structures is essential in many areas that are re- 
lated to high speed turbulent flow calculations. These include, for example, the 
predictions of high speed jet noise and chemically reacting flows. Therefore, this 
approach appears to be complementary rather than contradictory to the traditional 
second-order modeling techniques. 

6. Initiative to Develop New Second-Order Compressible Models 

To develop second-order models for compressible turbulence, a reasonable first 
step is probably to distinguish compressibility terms from incompressible back- 
ground. This is what Zeman (1990) and Sarkar et al. (1990) have performed. 
The dilatation dissipation they identified appears to be an important mechanism in 
compressible turbulence. A successful development of a compressible second-order 
model may inevitably invoke techniques and methodologies that were proved useful 
in the modeling of incompressible turbulence. The techniques and methodologies, 
however, have to be tuned in to the experimental observations and the theoreti- 
cal constrains of compressible turbulence. For example, one of the consequences 
of compressibility effects is the finite speed of the propagation of information. In 
supersonic flows, modulation of flow properties occurs only within Mach cones of 
influence with acoustic time delay. This feature introduces additional scales that 
may have to be included in the modeling of compressible flows. 

It can also be observed from the turbulent equations that there will be ex- 
changes between the turbulent kinetic energy and thermoenergy through rather 
different routes than in incompressible flows. In compressible turbulent flows of 
high fluctuating Mach numbers, eddy shocklets may begin to appear. Lumley in- 
dicated (private communication) that once the eddy shocklets appear in any sig- 
nificant numbers, a quite different approach, “... one which centers on entropy 
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fluctuations, and which concentrates on the rather large changes that take place in 
the shocklets...” is needed. Corresponding terms in the equations may need to be 
identified and properly modeled. These unique features, among others, may have to 
be considered explicitly in the development of new compressible turbulence models. 

Compressible turbulence modeling is still in its infancy. This seems to be true 
both theoretically and computationally. Therefore, it is prudent to have a open 
mind in the pursuit of rational, accurate and physically sound models. 
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